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ABSTRACT 

Presently  ice  extent  forecast  models  such  as  the  U.S.  Navy  Polar  Ice  Prediction 
System  (PIPS)  neglect  or  treat  small-scale  thennodynamic  processes  and  entrainment 
unrealistically.  Incorporating  better  algorithms  that  include  more  complete  physics  of  the 
mixed  layer  dynamics  will  allow  for  improved  prediction  of  ice  thickness  and 
distribution,  open  water  boundaries,  polynyas,  and  deep-water  fonnation  in  the  polar 
seas. 

A  one-dimensional  mixed  layer  turbulent  kinetic  energy  (TKE)  budget  model 
based  on  Garwood’s  NPS  mixed  layer  model  for  deep  convection  (Garwood,  1991)  was 
written  in  MATLAB®.  The  model  consisted  of  a  system  of  ten  equations  derived  by 
vertically  integrating  the  budgets  for  heat,  momentum,  salinity,  and  turbulent  kinetic 
energy  between  the  sea-ice-air  interface  and  the  base  of  the  turbulent  mixed  layer. 

The  NPS  mixed  layer  model  was  tested  using  atmospheric  forcing  and  ocean 
profiles  collected  at  the  Surface  Heat  Budget  of  the  Arctic  Ocean  Experiment  (SHEBA) 
site,  where  wind  stress  and  forced  convection  predominates.  Sensitivity  studies  using 
ocean  profiles  of  the  Greenland  Sea  were  also  conducted  to  address  thermodynamics  and 
ocean  profiles,  where  surface  cooling  and  free  convection  predominates. 

Specific  findings  and  results  include:  the  role  of  unsteadiness,  the  responses  of 
feedback  processes  depending  on  the  mixed  layer  structure,  and  the  importance  of  the 
temporal  resolution  of  the  model  forcing  on  both  skill  and  sensitivity  of  the  output.  The 
role  of  unsteadiness  such  as  heat  storage  within  the  mixed  layer  has  a  large  impact  on  ice 
melting  or  formation.  Feedback  between  the  atmosphere  and  ice  is  responsive  and 
depends  not  only  on  atmospheric  forcing  but  the  underlying  ocean  structure.  Feedback  is 
negative  for  stable  ocean  profiles  such  as  the  Western  Arctic  but  can  become  positive, 
promoting  deep  convection,  given  weakly  stable  ocean  profiles  like  that  of  the  Greenland 
Sea.  The  importance  of  model  time  step  is  evaluated  by  comparing  output  from  time- 
averaged  prescribed  surface  forcing.  The  long-term  goal  of  this  work,  and  future  studies, 
is  to  extend  the  relationships  determined  from  the  small  1-D  scales  to  the  larger  3-D 
scales  suitable  for  improvements  to  current  ice  models. 
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I.  INTRODUCTION 


The  scientific  objective  of  this  work  is  to  model  how  the  state  of  the  ice  cover 
changes  in  response  to  forcing  from  the  atmosphere  and  the  ocean.  Presently  ice  extent 
forecast  models  such  as  the  U.S.  Navy  Polar  Ice  Prediction  System  (PIPS)  neglect  or  treat 
small-scale  thermodynamic  processes  and  entrainment  unrealistically.  Incorporating 
better  algorithms  that  include  more  complete  physics  of  the  mixed  layer  dynamics  will 
allow  for  improved  prediction  of  ice  thickness  and  distribution,  open  water  boundaries, 
polynyas,  and  deep-water  fonnation  in  the  polar  seas.  As  a  result,  the  Navy  would 
improve  mission  planning,  effectiveness  and  efficiency  for  all  operations  in  the  polar 
regions.  Additionally,  because  the  same  physical  processes  that  are  explored  in  this  work 
as  mechanisms  for  ice  formation  and  melting  also  play  a  significant  role  in  the  global 
climate,  the  long-term  importance  of  this  study  goes  beyond  tactical  forecast  capabilities 
and  mission  effectiveness  to  include  macro-scale  circulation  and  deep-water  convection. 

The  sea-ice-air  system  in  the  high  latitudes  is  strongly  coupled.  Heat  and  salinity 
fluxes  that  determine  when,  and  how  much  ice  will  freeze  or  melt  are  a  result  of  the  sea- 
ice-air  system  interaction  as  a  whole;  no  single  factor  is  most  important.  Subtle  changes 
in  the  density  structure  and  circulation  of  the  polar  oceans  can  have  significant  effects. 
For  that  reason  differences  between  the  Western  Arctic  and  the  Greenland  seas  will  be 
addressed  in  order  to  shed  light  on  thermodynamic  properties  and  structures  that  enhance 
thennohaline  circulation,  thus  influencing  the  vertical  mixing  of  energy,  mass  and 
momentum  in  the  polar  seas.  Model  output  of  both  ice  thickness  and  ice  extent  can  be 
improved  by  including  the  atmospheric  processes  together  with  the  more  complete 
physics  of  the  mixed  layer  dynamics  and  the  underlying  temperature  and  salinity 
structure  of  the  ocean. 

A  one-dimensional  mixed  layer  turbulent  kinetic  energy  (TKE)  budget  model 
based  on  Garwood’s  NPS  mixed  layer  model  for  deep  convection  (Garwood,  1991)  was 
written  in  MATLAB®.  It  consisted  of  a  system  of  ten  equations  derived  by  vertically 
integrating  the  budgets  for  heat,  momentum,  salinity,  and  turbulent  kinetic  energy 
between  the  sea-ice-air  interface  and  the  base  of  the  turbulent  mixed  layer. 
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The  NPS  mixed  layer  model  was  tested  using  atmospheric  forcing  and  ocean 
profiles  collected  at  the  Surface  Heat  Budget  of  the  Arctic  Ocean  Experiment  (SHEBA) 
site.  Sensitivity  studies  using  ocean  profiles  of  the  Greenland  Sea  from  observations  by 
Quadfasel  and  Ungewip  (1988)  will  shed  light  on  the  subtle  thermodynamics  and  ocean 
profiles  that  enhance  thermohaline  circulation.  Findings  and  results  as  well  as 
recommendations  for  further  study  will  be  addressed  to  achieved  the  ultimate  goal,  to 
extend  the  relationships  determined  in  this  work  from  small  1-D  scales  to  the  larger  3-D 
scales  suitable  for  climate  models. 
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II.  POLAR  SEAS 


This  section  examines  the  physical  oceanography  of  the  polar  seas.  The  focus 
will  be  on  the  unique  temperature  and  salinity  structure  (Figure  1)  and  physical  processes 
that  influence  the  formation  and  melting  of  sea  ice.  The  sea  ice  cover  is  unique  to  the 
polar  oceans  and  affects  the  physiography  in  a  number  of  ways:  The  temperature  of  the 
surface  water  remains  near  the  freezing  point  as  a  function  of  its  salinity.  The  brine 
rejection  due  to  ice  freezing  increases  the  density  of  the  surface  waters  and  enhances 
thennohaline  circulation.  Atmospheric  winds  transfer  momentum  to  the  mixed  layer 
through  the  sea  ice  cover.  The  drag  coefficient,  and  thereby  the  momentum  transfer,  can 
vary  greatly  depending  on  the  character  of  the  ice,  e.g.  smooth,  rough.  Lastly,  the  ice 
albedo  is  seasonally  variable.  This  affects  the  exchange  of  solar  heating  at  the  surface 
and  the  energy  available  to  melt  the  ice  cover. 


Figure  1.  Polar  seas  unique  temperature  and  salinity  structure. 
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A.  ARCTIC  OCEAN  AND  ITS  MARGINAL  SEAS  CIRCULATION 
1.  Physiography 

Connected  primarily  to  the  Atlantic  Ocean,  a  lesser-input  circulation  is  from  the 
Bering  Strait,  which  provides  a  shallow  (less  than  45m)  connection  to  the  Pacific  Ocean; 
the  Arctic  Ocean  is  an  ocean  basin  that  occupies  an  area  of  approximately  12  million 
square  kilometers.  An  important  feature  of  this  ocean  basin  is  the  presence  of  sea  ice. 
Sea  ice  covers  less  than  10%  of  the  world’s  oceans,  and  40%  of  this  sea  ice  occurs  within 
the  Arctic  Ocean  basin.  The  presence  of  sea  ice  and  the  annual  cycle  of  melting  and 
freezing  have  a  significant  effect  on  the  exchanges  of  mass,  momentum  and  energy  of  the 
Arctic  Ocean  and  its  marginal  seas  (Figure  2). 
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Figure  2.  Topography  and  bathymetry  of  the  Arctic  (ET0P05  data  set,  NOAA 

1988). 
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The  temperature  of  surface  waters  in  the  Arctic  Ocean  basin  varies  from  -1.5  to 
-1.9  C  (i.e.,  the  freezing  point  of  sea  water  dependent  on  the  surface  salinity).  As  these 
surface  waters  move  through  the  Canadian  Arctic  Archipelago  and  Fram  Strait,  the 
temperature  increases  gradually.  Here  the  polar  surface  waters  exhibit  temperatures 
ranging  from  3  to  -1.9  C  (Figure  3). 

Surface  water  salinity  in  the  Arctic  Ocean  and  the  marginal  seas  is  relatively  low 
compared  to  other  oceans  (Figure  4).  In  the  Arctic  Ocean  itself,  surface  salinity  varies 
between  30  and  33,  and  decreases  in  the  area  of  the  shelf  seas  in  the  Western  Arctic  to 
below  30.  Generally,  the  salinity  is  lower  during  summer  than  winter  due  to  input  of 
freshwater  from  river  discharge  and  ice  melt.  During  the  winter,  the  combination  of 
sensible  cooling  and  brine  rejection  from  the  ice  increases  the  salinity  of  the  surface 
water  inducing  thennohaline  circulation  (AMAP  Assessment  Report,  1998). 
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Figure  3.  Winter  and  summer  surface  water  temperatures  (C)  in  the  Arctic  Ocean 
and  adjacent  seas  (USSR  Ministry  of  Defense  1980). 
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Figure  4.  Winter  and  summer  surface  water  salinity  in  the  Arctic  Ocean  and 
adjacent  seas  (USSR  Ministry  of  Defense  1980). 
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2.  Circulation 

The  Arctic  Ocean  water  mass  can  simply  be  described  as  waters  comprised  of 
relatively  warm  currents  originating  from  the  Atlantic  and  Pacific  Oceans  with  cold 
outflow  currents.  Atlantic  water  enters  the  Arctic  Ocean  through  Fram  Strait  and  the 
Barents  Sea,  while  a  small  percentage  of  Pacific  water  enters  via  Bering  Strait.  Water 
leaves  the  Arctic  largely  via  Fram  Strait,  but  also  through  the  Canadian  Arctic 
Archipelago  (MacDonald  and  Bewers,  1996)  (Figure  5).  Most  of  the  water  in  the  Arctic 
Ocean  originates  from  the  Atlantic  Ocean  (79%).  The  inflow  through  the  Bering  Strait  is 
a  small  percentage  of  the  total  (19%),  but  the  relatively  fresh  surface  water  inflow 
contributes  greatly  to  the  stability  and  to  the  sea  ice  extent  and  distribution  due  to  its 
stratifying  effect. 

The  main  water  outflow  is  via  the  East  Greenland  Current  (75%),  the  remaining 
outflow  (25%)  exits  via  the  Canadian  straits  (25%)  (AMAP  Assessment  Report,  1998). 


Figure  5.  Surface  ocean  currents  in  the  Arctic. 
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B.  MIXED  LAYER  STRUCTURE 

Differences  in  circulation  and  bathymetry  between  the  Arctic  Ocean  and  its 
marginal  seas  result  in  very  different  mixed  layer  structures.  Of  particular  interest  here  is 
the  contrast  between  the  Beaufort  and  Chukchi  Seas,  located  in  the  Western  Arctic,  and 
the  Greenland  Sea,  a  marginal  sea  located  south  of  the  Fram  strait  (Figure  2). 

1.  Beaufort  and  Chukchi  Seas 

The  mixed  layer  structure  in  the  Western  Arctic  is  complex  due  to  its  topography 
of  slopes,  ridges  and  deep-sea  plateaus.  Here,  waters  from  the  Pacific  and  from  the 
Atlantic  meet  forming  two  distinct  water  masses,  separated  by  a  cyclonic  boundary 
current  flowing  along  the  Northwind  Ridge  (Figure  6).  The  overall  impression  of  the 
Chukchi  and  Beaufort  Seas  is  an  ocean  with  considerable  horizontal  and  vertical 
structure,  and  with  a  variable  temperature  and  salinity  profile.  The  Western  Arctic  region 
is  highly  stable  due  to  the  inflow  and  mixing  of  the  buoyant  fresh  Pacific  Water  via  the 
Bering  Strait.  The  shallow  Chukchi  Sea  is  fed  primarily  by  Pacific  waters  entering  via 
the  Bering  Strait,  which  are  modified  and  mixed  with  Beaufort  Sea  water  in  the  interior 
of  the  Arctic  Ocean  along  the  Northwind  Ridge,  via  canyons.  The  Barrow  Canyon  is  a 
prime  location  for  such  interior  mixing.  The  boundary  current  on  the  Northwind  Ridge 
also  acts  as  a  source  of  interior  Arctic  waters  for  the  Chukchi  Sea  during  upwelling 
events  near  the  shelf  break  (Woodgate  et  al,  2004).  The  Pacific  inflow  to  the  Arctic 
Ocean  strongly  contributes  to  stabilizing  the  upper  ocean  in  the  Western  Arctic,  thereby 
influencing  ice  thickness  and  upper-ocean  mixing  (Aagaard  and  Carmack,  1994).  The 
inflowing  Pacific  water  is  less  dense  than  the  Arctic  Ocean  water  due  to  its  low  salinity 
and  higher  temperatures  compared  to  the  cold  Arctic  Surface  water.  This  creates  a  very 
fresh  mixed  layer,  over  a  saline  layer.  Because  salinity  dominates  the  density  profile  the 
temperature  excursions  are  less  constrained  by  buoyancy  forces  (Figure  7).  The  impact 
of  this  weak  and  variable  temperature  profile  is  that  thermohaline  circulation  does  not 
provide  a  significant  heat  flux  to  the  bottom  of  the  ice  in  the  Western  Arctic  (Woodgate 
et  al,  2004). 
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Figure  6.  A  schematic  of  the  circulation  over  the  Chukchi  Sea  and  Beaufort/Chukchi 

slope  (Weingartner,  2004) 
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T/S  profiles  from  SHEBA  CTD  Data 
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Figure  7.  Beaufort  and  Chukchi  Sea  variable  T/S  profiles  over  the  SHEBA  study. 

2.  Greenland  Sea 

The  circulation  between  the  Atlantic  and  the  Arctic  Ocean  results  in  a  very 
different  temperature  and  salinity  profile  in  the  Greenland  Sea  (Figure  8).  The  inflowing 
Atlantic  water  sinks  as  it  enters  the  Arctic  Ocean  basin  due  to  its  high  salinity  compared 
to  the  cold  but  fresher  Arctic  Surface  water.  This  Atlantic  and  Arctic  water  interaction 
creates  a  weaker  salinity  stratification  and  a  stronger  temperature  inversion  with  depth 
compared  to  the  Western  Arctic  (Figure  9).  Because  of  this,  the  Greenland  Sea 
thermohaline  circulation  readily  responds  to  cooling  and  ice  formation.  It  thus  entrains  a 
significant  amount  of  heat  from  below  and  transfers  heat  to  the  atmosphere  through  leads, 
polynyas  and  seasonally  open  seas,  allowing  large  surface  heat  fluxes.  When  ice  is 
present,  all  of  the  heat  entrained  into  the  bottom  of  the  mixed  layer  is  fluxed  upward  to 
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melt  the  ice  (Figure  10).  Therefore  the  sea-ice-air  exchanges  coupled  with  the  distinct 
ocean  structure  of  the  Greenland  Sea  promotes  deep-water  formation. 


Figure  8.  Representative  Greenland  Sea  T-S  profiles,  from  observations  by 

Quadfasel  and  Ungewip  (1988) 
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Western  Arctic  T/S  profile  comparison,  from  SFIEBA  (1997) 
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Figure  10. 
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When  ice  is  present  all  entrained  heat  is  fluxed  into  the  ice  (Garwood, 

2004). 
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C.  DEEP  CONVECTION 

Deep  convection  in  the  polar  seas  plays  a  key  role  in  the  Earth’s  overall  energy 
balance  by  exchanging  heat  stored  in  the  ocean  to  the  atmosphere.  It  is  understood  that 
strong  atmospheric  forcing,  wind  stress  and  surface  cooling,  must  be  present  in 
conjunction  with  other  conditions  for  deep  convection  to  occur  (Krauss  and  Businger, 
1994).  Strong  winter  cooling  of  northern  surface  waters  along  with  brine  rejection,  as  a 
result  of  ice  formation,  creates  dense  water  at  the  surface  that  sinks  and  mixes  with  water 
from  below.  It  causes  heat  to  be  transferred  from  the  warmer  ocean  to  the  colder 
atmosphere  or  melt  ice  when  present.  Ideally,  surface  water  must  have  a  relatively  high 
surface  salinity  so  that  the  surface  water  will  not  freeze  as  it  undergoes  the  intense 
cooling  necessary  in  order  to  become  dense  enough  to  sink  deeply.  The  integrated  heat 
and  salinity  fluxes  are  vitally  important  components  to  this  process  because  they 
determine  the  stability  of  the  upper  ocean.  The  short  time  scale  wind  stress  and  buoyancy 
fluxes  are  also  important  because  they  generate  turbulent  kinetic  energy  (TKE)  that 
causes  mixing  and  entrainment  thereby  deepening  the  mixed  layer  (Bramson,  1997). 

In  order  to  gain  a  better  understanding  of  the  important  processes  controlling  deep 
convection,  a  comparison  of  a  Western  Arctic  temperature  and  salinity  profile  and  a 
Greenland  Sea  wintertime  temperature  and  salinity  profile  will  be  used  in  conjunction 
with  representative  Arctic  atmospheric  forcing  from  SHEBA  to  illustrate  the  effects  that 
varying  temperature  and  salinity  profiles  have  on  deep  convection  under  the  same 
atmospheric  forcing  conditions. 

1.  Atmospheric  Forcing 

The  pattern  of  mean  sea-level  pressure  for  January  in  the  Arctic  shows  a  low- 
pressure  area  over  the  North  Atlantic  Ocean  around  southern  Greenland  and  Iceland 
(Icelandic  Low)  and  a  low-pressure  area  over  the  Pacific  Ocean  south  of  the  Aleutians 
(Aleutian  Low)  (Figure  11).  The  influence  of  the  Icelandic  Low  extends  through  the 
Fram  Strait  to  the  North  Pole.  Because  of  the  Icelandic  Low,  the  strongest  winds  occur  in 
the  Atlantic  sector  of  the  Arctic  where  they  follow  a  track  from  Iceland  to  the  Barents 
Sea.  The  Aleutian  Low  in  the  Pacific  is  blocked  by  the  mountains  of  Alaska  and 
northeast  Siberia  (Barry  and  Hare,  1974)  and  are  weaker  due  to  the  influence  of  the  high- 
pressure  areas  over  central  Canada  and  Tibet  (AMAP  Assessment  Report,  1998). 

15 


Figure  1 1 .  Mean  sea-level  pressure  for  January  in  the  Arctic  (AMAP  Assessment 

Report,  1998) 

In  July,  the  winds  weaken  throughout  the  Arctic  as  the  Aleutian  Low  disappears 
and  the  low-pressure  area  off  Iceland  shifts  to  southern  Baffin  Island  in  Canada 
(Canadian  Low)  (Figure  12)  (AMAP  Assessment  Report,  1998). 
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Figure  12.  Mean  sea- level  pressure  for  July  in  the  Arctic  (AMAP  Assessment  Report, 

1998). 

D.  NEAR  SURFACE  CONDITIONS  AND  OVERALL  ENERGY  BUDGET 

Near-surface  atmospheric  pressure  gradients  and  the  attributed  wind  fields  are 
primarily  responsible  for  forcing  the  ice  and  ocean  mixed  layer  circulation.  As 
mentioned  earlier,  this  stirring  entrains  heat  the  fluxes  into  the  underside  of  the  ice, 
resulting  in  melting.  If  the  vertical  circulation  is  weak  compared  to  the  stratification  of 
the  ocean  then  no  heat  will  be  entrained.  Additionally,  net  solar  radiation  can  act  to 
induce  melting  or  freezing.  Melting  ice  introduces  fresh  water  to  the  ocean  mixed  layer, 
creating  negative  salinity  flux,  and  inhibiting  entrainment.  Conversely,  ice  formation 
increases  the  salinity  at  the  surface  of  the  mixed  layer  due  to  brine  rejection,  positive 
salinity  flux,  and  enhances  entrainment.  The  formation  or  melting  of  ice  coupled  with  the 
mechanical  mixing  imparted  by  the  wind  stress  along  with  the  ocean  density  profile  will 
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determine  if  and  how  much  heat  will  be  entrained  and  provided  to  the  underside  of  the 
ice.  Ultimately,  ice  will  melt  or  freeze  as  a  result  of  the  overall  energy  balance. 

1.  SHEBA  Atmospheric  and  Ice  Measurement  Data 

Until  recently,  the  polar  region’s  atmospheric,  and  oceanographic  processes  have 
been  poorly  understood  due  to  lack  of  year-round  observations.  Due  to  the  short  time- 
scale  of  the  physical  processes  associated  with  sea-ice-air  interactions,  standard  climate 
observations  are  inadequate  for  mixed  layer  modeling  purposes.  Measurements  made  at 
SHEBA  provided  the  time-series  data  necessary  for  analyzing  and  understanding  the 
synoptically  and  seasonally  varying  atmospheric  and  oceanic  physical  processes. 

From  October  1997  to  October  1998,  the  SHEBA  ice  floe  drifted  more  than  1400 
km  in  the  Beaufort  and  Chukchi  Seas,  with  the  latitude  varying  from  74°N  to  81°N 
(Figure  13).  The  site  was  located  in  multiyear  pack  ice  with  summertime  meltponds  and 
occasional  nearby  leads. 


Figure  13.  SHEBA  ice  floe  drifted  more  than  1400  km  in  the  Beaufort  and  Chukchi 

Seas  (SHEBA  Map  of  Ice  Station  Drift,  2004). 
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The  drifting  ice  camp  scientists  (SHEBA  Principal  Investigators,  2004)  made 
measurements  of  the  atmospheric  boundary  layer,  ice  thickness  and  temperature  and 
salinity  profiles  of  the  underlying  ocean  (Figure  14-16)  (SHEBA  Phase  2  Data  Sets, 
2004). 


wind  speed  vs.  Time 


Figure  14.  SHEBA  atmospheric  data 
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Figure  16.  SHEBA  temperature  and  salinity  profiles  of  the  ocean  to  150  m. 

a.  Winds 

The  prevalent  wind  direction  was  from  40°  to  120°,  with  the  strongest 
winds  between  800-900  cm/s.  The  strongest  winds  occurred  in  the  winter  and  weakened 
during  the  summer  months.  Apparent  seasonal  changes  in  the  winds  recorded  at  the 
SHEBA  site  were  probably  due  to  a  combination  of  the  true  seasonal  wind  signal  and  the 
significant  change  in  geographic  location  due  to  the  ice  camp  drift  (Persson  et  al,  2002). 

b.  Total  Energy  Flux 

The  total  energy  flux,  Flot ,  into  the  surface  is  given  by 
Flol=Q0-Hs-Hl+C  (1) 
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where  C  is  the  conductive  flux,  Hs  is  the  sensible  heat  flux,  H,  is  the 
latent  heat  flux  and  Q0  is  the  total  net  radiative  flux  given  by 

Q0  =  -Qst  +  Qso  ~  Qu  +  Qio  (2) 

where  Qsi  is  incoming  shortwave  radiation,  Qn  is  incoming  longwave 
radiation,  Qso  is  outgoing  shortwave  radiation  and  Qlo  is  outgoing  longwave  radiation. 

The  assumption  in  (2)  is  that  all  radiative  flux  is  absorbed  within  the  ice.  This  is  a  good 
assumption  considering  our  ice  thickness  is  200-400cm  thick  (Persson  et  ah,  2002)  . 

For  simplicity,  the  surface  is  considered  as  a  slab  of  snow  and  ice  (Figure 
17).  Therefore  the  conductive  flux  (C)  is  implicit  in  the  radiative  flux  and  Ftot  becomes: 

Flol=Q0-Hs-H,  (3) 


turbulent 
mixed  layer 


Figure  17.  The  surface  is  considered  as  a  slab  of  snow  and  ice 
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Sensible  heat  flux,  Hs  and  latent  heat  flux,  II , ,  though  small  compared 
to  the  total  net  radiative  flux  Q0  (Figure  18),  are  important  because  the  sea-ice-air  system 
is  highly  sensitive  to  small  changes  in  the  net  heat  flux  due  to  the  coupled  nature  of  the 
system.  The  resulting  sign  convection  is  if  Ftot  is  negative,  the  heat  flux  is  into  the  ice. 

Thus  the  ice  is  gaining  energy,  which  acts  to  melt  the  ice.  If  Ftot  is  positive  this 
represents  cooling  and  more  ice  will  fonn. 

The  surface  energy  budget  for  SHEBA  shows  large  day-to-day  variability 
due  to  the  large  variability  in  the  sensible  and  latent  heat  flux,  shortwave  and  longwave 
radiation.  In  this  case  variations  in  the  net  solar  radiation  were  due  to  cloud  cover  effects 
and  albedo  fluctuations  due  to  changeable  surface  characteristics,  i.e.  fresh  snow  cover, 
white  ice,  and  meltponds.  The  net  shortwave  radiation  is  negative  from  March  to 
September,  and  the  net  longwave  radiation  is  positive  throughout  the  year.  The 
magnitudes  of  sensible  heat  flux  and  latent  heat  flux  are  a  factor  of  five  to  ten  times 
smaller  than  the  turbulent  heat  flux;  warming  the  surface  during  the  winter  and  weakly 
cooling  during  the  summer  time.  Variability  of  sensible  and  latent  heat  flux  is  owed  to 
the  atmosphere  vs.  surface  temperature  gradients,  which  is  affected  by  ocean  feedback 
due  to  the  presences  of  open  water  leads,  and  wind  speed  (Figure  19).  Summing  the 
radiation  flux,  sensible  and  latent  heat  flux  yields  a  net  heating  from  May  through  August 
and  cooling  for  the  period  of  September  through  March  (Persson  et  ah,  2002). 
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Figure  18.  Heat  Flux  Components 
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Figure  19.  Sensible  heat  flux  and  latent  heat  flux  warming  the  surface  during  the 
winter  and  weakly  cooling  during  the  summer  time 

c.  Ice  Thickness 

The  spatial  and  temporal  evolution  of  sea  ice  is  a  combination  of  dynamic 
and  thermodynamic  processes.  The  observed  thickness  distribution  is  a  result  of  the 
coupled  mechanism  of  thermodynamic  ice  growth  or  melting  and  the  compaction  and 
dilation  of  the  ice  due  to  the  dynamics  of  wind  stress  and  ocean  currents,  which  modify 
the  spatial  distribution  of  sea  ice  (Flarder  and  Lemke,  1994).  Daily  estimates  of  the  ice 
velocity,  deformation,  mean  thickness,  and  thickness  distribution  were  taken  during  the 
SHEBA  experiment.  From  the  start  of  the  experiment  in  October  to  early  mid  June,  165 
cm  of  ice  fonned.  The  summer  melt  season,  lasting  52  days  from  June  through  August, 
resulted  in  108  cm  of  ice  melted.  After  this,  ice  formation  continued;  17  cm  was 
observed  from  early  August  to  the  conclusion  of  the  experiment  in  October  1998.  The 
data  shows  a  net  gain  of  74  cm  over  a  one -year  period.  This  is  not  taken  to  be  a  result  of 
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thermodynamic  processes.  It  can  probably  be  attributed  to  the  change  in  geographic 
location  due  to  the  ice  camp  drift  and  dynamic  processes  such  as  convergence  of  the  ice. 
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III.  PROCEDURE 


A.  MIXED  LAYER  ICE  COUPLE  MODEL  (MICE) 

1.  Ocean  Mixed  Layer  Structure 

The  idealized  ocean  mixed  layer  is  a  homogeneous  fully  turbulent  region  of  the 
upper  ocean  that  is  bounded  above  by  the  air-sea  interface,  and  below  by  a  turbulent 
entrainment  zone.  This  entrainment  zone  is  where  temperature  and  salinity  undergo 
dynamic  instabilities,  which  give  rise  to  entrainment  fluxes  of  temperature  and  salinity. 
Below  this  entrainment  zone  is  stratified  water  that  increases  in  density  with  depth.  The 
source  of  energy  for  the  generation  of  turbulence  in  the  mixed  layer  is  provided  by  the 
fluxes  of  heat,  salt  and  momentum  at  the  sea-ice-air  interface  (Garwood,  1977). 

2.  Turbulent  Kinetic  Energy  Budget 
a.  Notation 

Seawater  variables  that  vary  in  time  and  space  may  be  decomposed  into 
mean  and  fluctuating  parts.  For  example,  total  potential  temperature  may  be  expressed 
as 0(x,y,z,t)  =  0(z,t)  +  6'(x,y,z,t) ,  where  the  overbar  is  a  horizontal  mean  over  a 
distance  or  horizontal  model  grid  size  L  that  is  large  compared  with  the  scale  of  the 


fluctuations  (turbulence  integral  scale).  8  -  T .  All  mixed  layer  variables  may  be 


decomposed  similarly  into  mean  and  turbulent  parts,  T  =  T  +  T ' , 
S  =  S  +  S'  ,u  =  w  +  w',v  =  v  +  v',  and  w  =  w  +  w' . 


In  addition  to  the  horizontal  mean,  there  is  a  vertical  mean  denoted  by 


—  1  '  — 

that  is  a  vertical  average  over  the  mixed  layer  of  depth  h  ,  l6  ''j  =  —  J  ddz  between  the 


bottom  of  the  mixed  layer,  z  =  —h  and  the  bottom  of  the  ice,  z  =  -ht . 


Turbulent  fluxes  are  associated  with  velocity-property  covariances.  For 
example,  vertical  heat  flux  is  pCp6'W'  (Garwood,  2004). 
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b.  Entrainment  Hypothesis 

The  TKE  budget  is  the  basis  for  the  entrainment  hypothesis.  The 
assumption  is  that  turbulence  contained  within  the  mixed  layer  provides  the  energy 
needed  to  destabilize  and  erode  the  underlying  stratified  water  mass.  The  amount  of  heat 
provided  to  the  underside  of  the  ice  is  dependent  on  this  understanding  of  the  dynamics 
of  the  mixed  layer  deepening  or  retreat  (Garwood,  1977).  The  one-dimensional  total 
TKE  equation  is: 
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where  E  =  u'2  +v'2  +w'2 ,  u  is  the  eastward  velocity  component,  v  is  the 
northward  velocity  component,  and  w  is  the  upward  velocity  component  (z  is  positive 
up).  To  simplify  the  system  of  equations,  horizontal  homogeneity  is  assumed.  The 
approximation  of  local  homogeneity  holds  because  of  the  short  time  scale  for  turbulence. 
Garwood  (1977)  recognized  that  the  mean  buoyancy  and  momentum  fluxes  are  not  one¬ 
dimensional  for  all  time  and  space  scales  but  effects  such  as  advection  are  neglected  in 
order  to  focus  on  new  aspects  of  mixed  layer  modeling  such  a  buoyancy  flux 
enhancement.  The  terms  of  (4)  are:  time  rate  of  change,  shear  production  of  horizontal 
TKE,  buoyancy  flux,  transport  of  turbulence  and  viscous  dissipation. 
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:  Time  rate  of  change  or  storage  terms.  In 


steady  state  these  terms  are  equal  to  zero.  These  terms  can  be  positive  or  negative  but 
due  to  the  rapid  adjustment  of  turbulence  they  are  assumed  to  be  negligible  compared  to 
the  production  and  dissipation  terms. 


dz 
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Shear  production  of  horizontal  TKE. 


The 


production  of  forced  convection  due  to  wind  stirring.  This  term  can  only  be  positive  and 
is  a  large  source  of  TKE. 


b  V  :  Buoyancy  Flux.  This  term  can  be  positive  or  negative.  When 
buoyancy  flux  is  positive  it  is  a  source  of  vertical  TKE.  An  unstable  density  profile  will 
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result  in  free  convection  and  positive  buoyancy  fluxes.  However,  if  the  oceanic  profile  is 
stable  the  buoyancy  flux  will  damp  the  mechanical  effects  of  the  shear  production  and 
convert  TKE  to  potential  energy  (PE). 

:  Transport  of  turbulence.  Turbulence  is 

advected  when  w'  carries  TKE  or  when  fluid  parcels  collide  passing  one  parcel’s 
momentum  and  kinetic  energy  (KE)  to  the  other.  The  advection  and  pressure  transport 
tenns  do  not  change  the  total  energy  in  the  mixed  layer  they  only  move  it  around. 


d_ 

dz 


f  u'2+v'2+w'2  p 
w'\ - +  — 


V 


PoJ 


s  —  3s  -  3 £  -  3 s  ;  Viscous  Dissipation.  This  is  a  destruction  term 

and  can  only  be  negative.  At  the  small  time  and  space  scales  the  turbulence  is  isotropic. 
Therefore  each  of  the  different  dissipation  components  is  equal  to  each  other  and  equal  to 
one  third  of  the  total  dissipation  (Garwood,  1977). 
c.  Integration  of  TKE 

Assuming  steady  state,  and  horizontal  homogeneity  the  TKE  budget  (4) 
may  be  integrated  from  the  bottom  of  the  mixed  layer,  z  =  —h  to  the  bottom  of  the 
ice,  z  =  -ht . 


dl^jh 

dt 


1 


■dU 

dz 


■dV 

dz 


(5) 


The  change  in  TKE  in  the  mixed  layer  can  be  related  to  entrainment  in 
tenns  of  stratification  and  mixed  layer  depth,  h.  Stratification  relates  entrainment 
velocity  and  mean  property  jumps  between  the  bottom  of  the  mixed  layer  and  the 
underlying  layer  to  the  vertical  flux  of  the  property  out  of  the  entrainment  zone  into  the 
base  of  the  mixed  layer.  Entrainment  velocity  can  be  described  by  the  equation  of  mixed 
layer  depth: 


dh  — 
—  =  w-W 
dt 


(6) 


Assuming  no  vertical  motion,  W 


h 


=  0 ,  the  equation  becomes: 
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(7) 


dh 


Wind  stress  at  the  surface  produces  turbulence  in  the  surface  layer  this 
turbulence  deepens  the  mixed  layer  mixing  temperature,  salinity,  and  momentum  into  the 
mixed  layer. 


=  -weA6  =  heat  flux 


S'W 


=  -weAS  =  salinity  flux 


U'W1 


=  -w  AU  =  momentum  flux  of  eastward  momentum. 


From  this  forced  convection  can  be  solved  as  a  combination  of  the  shear 
production  at  the  surface  and  entrained  momentum  flux  at  the  bottom  of  the  mixed  layer. 


-h  L 


■dU 

dz 


u'w' - l -vV 


8V 

dz 


dz  =  2m3ul  +  ((u^J  +  (v^  )w( 


(8) 


where  ul  is  a  wind  forcing  scale  and  m3  is  a  shear  production  coefficient. 


The  buoyancy  flux  term  can  be  expressed  in  terms  of  heat  and  salinity 
fluxes  by  the  use  of  the  linearized  equation  of  state: 


p  =  po[l-a(0-0o)  +  fi(S-So)] 


(9) 


The  buoyancy  is  defined  as: 

b7w'  =  ag07w'-j3g'S7w'  (10) 

where  a  is  the  thermal  expansion  coefficient  and  ft  is  the  salinity 
contraction  coefficient. 

Due  to  the  pressure  dependence  of  the  thermal  expansion  coefficient 
buoyancy  is  not  conserved.  Unlike  forced  convection,  buoyancy  flux  cannot  be 
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expressed  as  a  simple  combination  of  the  surface  buoyancy  flux  and  mixed  layer 
buoyancy  flux.  The  vertical  integral  of  the  buoyancy  flux  must  be  calculated  (Garwood 
1991).  The  thermal  expansion  coefficient  varies  with  temperature,  salinity  and  pressure, 
the  thermobaric  effect.  This  effect  is  small  for  shallow  mixing  but  can  be  significant  in 
deep  mixing  cases.  The  variation  with  depth  is  nearly  linear  (Figure  20). 


Representative  Western  Arctic  Alpha  vs.  Pressure 


Representative  Greenland  Sea  Alpha  vs.  Pressure 


Figure  20.  Linear  variation  of  the  thermal  expansion  coefficient  with  temperature, 

salinity  and  pressure 

The  thermal  expansion  coefficient  can  be  written  as: 

a~ao(0,S)-ax(0,S)z  (11) 

where  a0  is  the  thermal  expansion  coefficient  at  the  surface,  and  ax  is  a 

coefficient  for  the  pressure  enhancement  of  buoyancy.  To  add  physical  meaning  to  the 
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importance  of  the  pressure  enhancement  of  buoyancy  flux  Garwood  (1991)  defined  a 
“thermobaric  depth”  at  which  the  increased  pressure  doubles  the  thennal  expansion 
coefficient  relative  to  its  surface  value. 


(12) 


The  thennobaric  depth  increases  strongly  with  increasing  potential 
temperature  and  weakly  with  increasing  salinity  (Figure  21). 


Temp 


Figure  21.  Dependence  of  Ha  on  potential  temperature  and  salinity 
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Garwood  (1991)  defines  a  dimensionless  depth  of  mixing  as  - . 

H a 

h 

Buoyancy  is  conserved  for  - «  1 .  In  warmer  ocean  waters  the  thennobaric  effect  is 

Ha 

negligible  to  depths  of  about  1000m.  However,  in  cold  polar  oceans  the  dependence  of 
a  upon  pressure  cannot  be  neglected  for  mixing  depths  greater  than  100m  because  of 
potential  temperature’s  effect  on  the  thennobarric  depth  (Garwood  1991). 

In  terms  of  thennobaric  depth  (12)  the  thermal  expansion  coefficient 

becomes: 


f 

a  =  ag  1 

V 


(13) 


Buoyancy  flux  combined  with  equation  (13)  yields 


b'w'  =  b'w'g  —b'w's 


(  7  ^ 

“ 

( 

b'w'e  =  ao 

{\~  — 

l  Ha) 

gO'w'  =  ag 

0'w'  0 

1  +  - 

l  h ) 

-0'w'. 

(14) 


b'w\  =  j3gS'w'  =  J3ag 


( 

(  vW 

5'w'0| 

1  +  - 

l  h ) 

-S'w'-„\ 

i  ^ 

UJ. 

The  vertical  integral  of  buoyancy  flux  can  now  be  solved  for  deep  ocean 
mixing  and  polar  seas  cases. 


u  u 

J  b'w'^z^dz  =  J  [b'wg '  +  b'ws  '^jdz 


Taking  the  vertical  integral  of  the  potential  temperature  component  of 
buoyancy  flux,  b'we'\ 
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Vertical  integral  of  the  salinity  component  of  buoyancy  flux,  b'w s 


a  -wr-,  ,  S'w'oh  S'w'-hh 
=  /3og  S'w'0h - — + — - - 


=  PoE 


S'w'oh  S'w'-hh 

- 1 - 

2  2 


The  resulting  equation  (15)  illustrates  the  effect  of  the  vertical  integral  of 
buoyancy  flux,  and  its  dependence  on  the  vertical  entrainment  velocity.  The  potential 
temperature  component  of  buoyancy  flux  is  dependent  on:  surface  temperature  flux, 
temperature  at  the  bottom  of  the  mixed  layer,  mixed  layer  depth,  and  thennobaric  depth. 
Since  the  haline  contraction  coefficient  is  not  a  function  of  pressure,  the  salinity 
component  of  buoyancy  flux  is  only  dependent  on:  surface  salinity  flux,  (precipitation 
and  evaporation  with  no  ice  present,  melting  and  freezing  when  ice  is  present),  salinity  at 
the  bottom  of  the  mixed  layer,  and  mixed  layer  depth. 
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(15) 


J  b'w'fz^dz  = 


fh  h2  'l 

fh  h2  ) 

~ 

( h > 

f  h\] 

G'w\ o 

+ 

+  0'w'-h 

+ 

+  jB0g 

S'w'o 

+  SV-i 

l2  6HaJ 

l2 

^21 

3.  1-D  Mixed  Layer  Ice  Coupled  Model 

The  Mixed  layer-ICE  coupled  model  (MICE)  is  a  one-dimensional  second-order 
turbulence  closure  solution  for  the  first  and  second  moments  of  the  unsteady  oceanic 
planetary  boundary  layer,  including  thermobaric  enhancement  of  entrainment  (Garwood, 
1991).  The  MICE  coupled  model  has  ten  prognostic  equations  obtained  by  integrating 
the  mixed  layer  over  the  depth. 

Mixed  Layer  Depth 


dh 


(16) 


Ice  Thickness 


—L  =  F  +  Pr-Ev  (17) 

dt 


When  ice  is  present,  precipitation,  Pr,  and  evaporation,  Ev,  are  zero.  F,  is 
the  formation  or  melting  of  ice: 


F=(Q  0-pCP6'W' 


y(Pi  l.) 


where  Q0  is  the 


flU-X,  Q0  Q Back  Q- Sensible  O' Latent 

Heat 


-'Back 
Radiation 


-'Latent 

Heat 


-'Solar 

Radiation 


and 


net 


upward  surface  heat 
is  the  entrainment  heat  flux  into 


the  bottom  of  the  ice  from  the  ocean. 
Integrated  total  TKE 


d 


E)h  _ . 

-^  =  G0  +  G/l+(B'W')h+D 


(18) 
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is  the  sum  of  the  vertically  integrated  shear  production  (8)  and  buoyancy  flux  (15),  plus 
viscous  dissipation,  s .  An  estimation  of  viscous  dissipation  can  be  made  by  taking  the 
rate  at  which  the  largest  eddies  supply  energy  to  the  smaller  eddies  to  be  proportional  to 
the  reciprocal  of  the  time  scale  of  the  largest  eddies.  The  net  rate  of  dissipation  and  the 
mean  turbulent  energy  are  used  to  define  the  dissipation  time  scale,  re  (Garwood,  1977). 


If  the  time  scale  of  the  largest  eddy  is  proportional  to  the  mixed  layer  depth 
divided  by  the  rms  turbulent  velocity, 


h 


then, 


o  _ 

D  =  J  <?dz  =  -2m,  (e 

-h 

where  mi  is  the  dissipation  coefficient. 
Integrated  vertical  TKE  equation  is: 


3/2 


W,2)h 


-^  =  {B,W'jh+Rp+D/3 


(19) 


Rp  is  the  redistribution  of  turbulent  energy  defined  as: 
-3^W'2|)j^E^  (Lumley,  1974).  The  rational  is  that  of  a  “return  to 

isotropy”,  meaning  the  pressure/strain  interaction  will  tend  to  distribute  equally  energy 
among  u,  v,  and  w  components. 

The  vertical  integral  of  potential  temperature  flux  from  the  base  of  the  mixed 
layer  to  the  sea-ice  interface  defines  the  change  in  mixed  layer  potential  temperature 


where 
RP=  2m2 
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A  n\  -O'W'  -O'W' 

\  -h,  -h 


Recall,  O'W'  =  -weAO  .  Thus  the  change  in  mixed  layer  salinity  can  be  defined 


in  the  same  way,  or  S'w  =  -weAS  . 


d Is)  -s'ff"  .  -sw 


The  temperature  flux  at  the  top  of  the  mixed  layer 


O'W'  =CTCDllNE[(0)-Tfp 


must  be  physical  linked  to  the  freezing  point  of  the  surface  water,  which  is  a  function  of 
the  surface  salinity  flux,  i.e.  melting  or  freezing  of  ice. 


1,-#K 


This  interaction  between  temperature  flux  and  salinity  flux  illustrates  the  strong 
coupling  within  the  sea-ice-air  system. 

Mixed  layer  current  and  ice  speed  are  calculated 


J^  =  2Qz(F)^  +  ^-(f/)we-Y^ 


dt  z\  /  p  \  /  Tm 


dLlA  =  2nzvk+—-—-— 

dt  '  Pi  Pi  Tm 


V  h  r--  V  h 

d  _rh_  _  _2Q  U ,h,  +  ^  ^  -  M. 

dt  Pi  Pi  Tm 
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where  Qz  =  7.29x10  5  sin(/at)  is  the  Coriolis  force.  tLx  is  the  atmospheric  wind  stress 
acting  on  the  ice  surface  in  the  x-direction.  r0v  is  the  stress  at  the  sea-ice  interface  acting 
on  the  ocean  surface  in  the  y-direction.  Tm  is  a  damping  term  to  dissipate  energy  within 
the  system. 

The  entrainment  velocity  is  calculated  by  vertically  integrating  the  TKE  budget 
(Bramson,  1997),  giving: 

we  =  (w'2  y/2 

This  gives  a  system  of  ten  prognostic  equations:  mixed  layer  depth,  (26),  ice 
thickness,  (27),  total  TKE,  (28),  vertical  TKE,  (29),  mixed  layer  potential 
temperature,  (30),  mixed  layer  salinity,  (31),  mixed  layer  velocity  components,  (32) 
and(33),  and  ice  velocity  components,  (34)  and  (35);  for  a  sea-ice-air  mixed  layer  model 
that  can  be  solved. 
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IV.  SENSITIVITY  STUDIES 


The  following  sensitivity  studies  demonstrate  the  strongly  coupled  feedback 
system  for  sea-ice-air  interactions  in  the  high  latitudes.  Realistic  high  arctic  synoptic 
forcing,  provided  by  SHEBA  wind  and  net  radiation  fluxes,  were  used  to  examine  the 
effects  of  MICE  model’s  sensitivity  to  varying  initial  conditions.  Initial  conditions 
include:  turbulent  kinetic  energy,  mixed  layer  depth,  and  temperature  and  salinity 
profdes. 

Representative  Western  Arctic  and  Greenland  Sea  profiles  were  used  to 
demonstrate  how  important  realistic  ocean  structure  is  to  forced  convection,  free 
convection  and  total  TKE.  The  dependence  of  ice  melting  or  formation  on  the  magnitude 
of  wind  speed  was  studied  by  varying  wind  forcing  by  50%  to  100%  of  the  observed 
forcing  with  the  same  surface  cooling  and  same  Greenland  Sea  temperature  and  salinity 
profile.  The  effect  of  the  net  heat  flux  was  studied  by  running  simulations  with  heat  flux 
values  that  were  about  +/-  0.0002  cal  cnf“  s'  of  the  observed  forcing,  with  the  same  wind 
forcing  and  same  Greenland  Sea  temperature  and  salinity  profiles.  The  decrease  or 
increase  in  net  heat  flux  was  representative  of  the  warming  attributable  to  the  increased 
sensible  heating  due  the  presence  of  warmer  air  aloft  or  cooling  effects  attributable  to 
latent  heat  cooling  due  to  the  presence  of  leads.  Finally,  a  moving  average  was  applied  to 
the  atmospheric  forcing  to  simulate  daily,  weekly  and  monthly  forcing.  The  simulated 
daily,  weekly,  and  monthly  atmospheric  data  along  with  the  measured  hourly  data  was 
used  to  examine  the  effects  on  the  MICE  model’s  skill  due  to  the  varying  time  scales  of 
the  forcing  data. 

The  values  for  the  physical  constants  are  given  in  Table  1  and  surface  boundary, 
and  below-layer  forcing  conditions  for  the  Western  Arctic  and  Greenland  Sea  are  given 
in  Table  2. 
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Description 

Symbol 

C-G-S  Value 

Freezing  Temperature 

T 

Lfp 

Function  of  Salinity 

Surface  thermal  expansion  coeff. 

°o 

7.44x1  O^C'1 

Thermobaric  depth 

Ha 

215220  cm 

Haline  contraction 

00 

7.94x1  O'4 

Ice  salinity 

s, 

6  psu 

Representative  seawater  density 

P 

1.028  gm/cm3 

Ice  density 

Pi 

0.905  gm/cnr 

Representative  air  density 

Pa 

.00125  gm/cm3 

Latent  heat  of  fusion  for  seawater 

A 

66.88  caPgm 

Heat  capacity  of  seawater 

cP 

0.9297  caFgm/C 

Ice  thermal  conductivity 

A 

0.0048  cahs/cm/C 

Drag/heat  transfer  coefficients 

cD  =  cT 

1.5xl0'3 

Planetary  rotation 

Qz 

1.4x10V 

Dissipation  coeff. 

m] 

1 

Pressure  Redistribution,  coeff. 

m2 

0.5 

Shear  Production  coeff. 

m. 

6 

Damping  Coefficient 

T 

m 

8640  s 

Table  1.  Table  of  Values  for  Physical  Constants. 
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Description 

Symbol 

C-G-S  Value 

Surface  heat  flux 

9) 

SHEBA  Measured  Data 
(Cal/cnr/s) 

Eastward  wind  stress 

Tx 

Calculated  (dyn/cm  ) 

Northward  wind  stress 

Ty 

Calculated  (dyn/cm  ) 

Vertical  Mean  Velocity  at  bottom  of  mixed  layer 

W-h 

0  cm/s 

Western  Arctic  Temperature  just  below  mixed 
layer 

0  below 

-1.07  C 

Western  Arctic  Salinity  just  below  mixed  layer 

S  below 

29.31 

Greenland  Sea  Temperature  just  below  mixed 
layer 

0  below 

-0.33  C 

Greenland  Sea  Salinity  just  below  mixed  layer 

S  below 

34.88 

Table  2.  Table  of  Surface  Boundary  and  Below-Layer  Forcing  Conditions  . 


A.  VARYING  INITIAL  CONDITIONS 

1.  Turbulent  Kinetic  Energy 

Sensitivity  to  varying  initial  TKE  was  examined  by  running  the  MICE  Model 
with  purposefully  too-high,  reasonable,  and  too-low  initial  TKE  values.  For  the 
Greenland  Sea  case  values  of  0.1  cm  /s  ,  0.3  cm  /s'  and  1.0  cm'/s  ,  for  low,  best  guess 
and  high  values  respectively  were  used  (Figure  22-23).  Examining  the  resulting  plots  of 
TKE  shows  the  TKE  converging  quickly  to  a  balanced  solution.  The  model  achieves  a 
good  solution  within  10  hours,  and  the  TKE  solution  is  fully  spun  up  by  0.5  day  and  does 
not  diverge  throughout  the  model  run. 
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s/  wo 


Figure  22.  Sensitivity  to  varying  TKE  -  Greenland  Sea  Julian  Day  308-318. 
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Greenland  Sea  Varying  Initial  Energy  E  &  w2  vs.  Time 


Figure  23.  TKE  Convergence  at  0.5  Day-  Greenland  Sea  Julian  Day  308-318. 

2.  Mixed  Layer  Depth,  Temperature  and  Salinity  Profile 

Western  Arctic  and  Greenland  Sea  profiles  were  run  using  the  same 
representative  synoptic  forcing  to  illustrate  the  effect  of  the  initial  underlying  ocean 
structure  on  ice  melting  and/or  formation.  Both  cases  were  run  under  relatively  light 
wind  conditions  where  the  circulation  was  primarily  thermally  driven  due  to  surface 
cooling  and  ice  freezing.  The  Western  Arctic  initial  mixed  layer  depth  was  2800  cm, 
initial  mixed  layer  temperature  was  -1.51  C  and  initial  mixed  layer  salinity  was  27.75. 
The  Greenland  Sea  initial  mixed  layer  depth  was  199  m.  Initial  mixed  layer  temperature 
was  -1.91  C,  based  on  the  freezing  point  at  the  initial  mixed  layer  salinity  of  34.79. 
Figure  24  is  a  comparison  of  the  two  different  initial  profiles. 
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.^4  Representative  Western  Arctic  and  Greenland  Sea  Temperature  and  Salinity  Profiles 


Figure  24.  Representative  Ocean  Structure  for  Western  Arctic  and  Greenland  Sea. 

The  Greenland  Sea  has  weaker  salinity  stratification  and  a  stronger  temperature 
inversion  than  does  the  Western  Arctic.  Therefore,  by  the  entrainment  hypothesis,  we 
would  expect  more  mixing.  The  results  show  this  to  be  the  case  (Figure  25).  During  the 
fifty  day  MICE  model  run  the  representative  Greenland  Sea  mixed  layer  deepens  2000 
cm  whereas  the  representative  Western  Arctic  mixed  layer  deepens  500  cm. 

Due  to  the  net  surface  cooling  there  is  ice  freezing  in  both  cases,  but  the  greater 
entrained  heat  into  the  ice  from  the  underlying  ocean  in  the  Greenland  Sea  profile  results 
in  less  ice  formation.  The  Greenland  Sea  forms  5.5  cm  of  ice  during  the  fifty  day  run. 
Under  the  same  atmospheric  conditions  the  Western  Arctic  forms  35.5  cm.  The  salinity 
increase  of  0.51  is  associated  with  the  physical  process  of  brine  rejection  due  to  ice 
freezing.  As  the  salinity  increases,  the  freezing  point  decreases.  Figure  26  shows  heat 
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storage  for  both  the  Greenland  Sea  and  Western  Arctic  mixed  layer.  Physically,  the  ice  is 
acting  as  an  insulating  cap  that  does  not  allow  entrained  heat  to  directly  escape  to  the 
atmosphere.  Recalling  the  prognostic  equation  for  the  mixed  layer  temperature  is 


-O'W' 

-hj 
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A6*we 

- ,  the  MICE  model  allows  of  this  for  the  heat  storage  by 


regulating  the  heat  flux  into  the  ice.  The  temperature  flux  at  the  top  of  the  mixed  layer  is 
-,/2  /E((0)-W 
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fP 


where  (6)-Tfp  is  the  difference  of  mixed  layer 


temperature  from  the  freezing  point.  Fluctuations  in  this  temperature  difference  are  an 
order  of  magnitude  smaller  than  the  fluctuations  in  TKE.  The  square  root  of  Total  TKE 
present  in  the  system,  \J~E ,  governs  the  transfer  of  heat  into  the  ice  and  thereby  heat 
storage  in  the  mixed  layer  (Figure  27). 


Change  in  Mixed  Layer  Depth  vs.  Time 


Change  in  Ice  Thickness  vs.  Time 


Figure  25.  Sensitivity  to  Initial  conditions:  h,  T(z),  S(z). 
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Heat  Storage  (T-Tfp) 
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Figure  26.  Fleat  storage  for  both  the  Greenland  Sea  and  Western  Arctic  mixed  layer. 
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Total  TKE  vs.  Time 


T-Tfp 


Julian  Day 


Figure  27.  Difference  of  mixed  layer  temperature  from  the  freezing  point  is  an  order 

of  magnitude  smaller  than  the  TKE  in  the  system 

B.  VARYING  WIND  SPEED  WITH  THE  SAME  NET  HEAT  FLUX 

The  dependence  of  melting  and  freezing  on  wind  speed  was  examined  by  running 
three  different  wind  forcing  cases  with  the  Greenland  Sea  temperature  and  salinity 
profiles  (Figure  24)  for  thirty  days.  The  SHEBA  representative  synoptic  wind  timeseries 
(Figure  14)  was  used  for  case  1.  As  a  control,  the  SHEBA  winds  were  applied  to  the 
Greenland  Sea  case  in  order  to  isolate  the  dependence  upon  initial  ocean  conditions  using 
realistic  synoptically  varying  winds.  Case  2  the  same  wind  was  increased  50%  and  for 
case  3  the  wind  forcing  was  increased  100%. 

The  results  illustrate  the  complex  and  coupled  nature  of  the  sea-ice-air  system 
(Figure  28).  The  system  melts  and  freezes  ice  in  pulses,  with  significant  negative 
feedback  tied  in  with  the  periods  of  rapid  entrainment.  Generally  the  deepening  mixed 
layer  entrains  heat,  which  increases  the  mixed  layer  temperature  and  heat  storage.  Then 
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the  extra  heat  contained  in  the  mixed  layer  is  transferred  to  the  underside  of  the  ice  in 
proportion  to  the  square  root  of  total  TKE  present  in  the  system,  \f~E  .  At  the  same  time 
the  salinity  fluxes  due  to  ice  melting  and  fonnation  impact  the  buoyancy  flux  which  has  a 
direct  effect  on  the  entrainment. 


The  mixed  layer  depth  plot  (Figure  28)  shows  slightly  greater  deepening  for  the 
lowest  wind  case  and  the  least  deepening  for  the  highest  wind  speed  case.  This  result  is 
due  to  the  effect  on  free  convection  of  the  salinity  flux  from  ice  melting.  The  melting  ice 


creates  a  positive  salinity  flux,  SW' 
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at  the  top  of  the  mixed  layer,  which 


in  turn  damps  the  free  convection  by  inducing  a  negative  buoyancy  flux, 
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The  degree  to  which  buoyant  damping  or  enhancement  is  occurring  is  very 
complex  because  simultaneously  the  increasing  wind  is  providing  more  heat  to  the  ice 
bottom  due  to  the  increased  TKE  (Figure  29).  The  greater  heat  to  the  ice  bottom  is 
apparent  in  the  change  in  ice  thickness  between  the  three  cases.  Case  1  and  case  2  show 
a  modest  ice  formation  of  1.22  cm  and  0.55  cm  respectively.  However  the  highest  wind 
case,  case  3,  resulted  in  0.17  cm  of  ice  melting  despite  the  smallest  entrainment  values. 
Due  to  the  increased  TKE  of  the  highest  winds,  the  heat  storage  in  the  mixed  layer  was 
less,  and  more  heat  was  transferred  into  the  ice  bottom  from  the  underlying  ocean. 
Figure  30  shows  the  high  frequency  fluctuations  at  which  ice  melting  or  freezing  will 
occur.  Ice  freezing  is  based  on  the  net  atmosphere  and  ocean  inputs.  Heat  is  lost  to  or 
gained,  from  the  atmosphere  and  heat  lost  to  or  gained,  from  the  ocean.  Ice  losing  heat  to 
the  ocean  seems  unphysical  but  if  the  freezing  point  of  seawater  is  less  than  the 
temperature  of  the  ice,  the  ice  can  give  up  heat  to  the  water. 


F=(Q  0-PCP0'W' 


)/(Pi  Ls) 


If  heat  is  lost,  more  ice  will  form.  If  heat  is  gained,  the  ice 


will  melt. 
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Figure  28.  Varying  Wind  Speed  same  Net  Heat  Flux 
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Total  TKE  vs.  Time 
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Figure  29. 


Total  TKE  and  vertical  TKE  for  vary  wind  speed. 
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Figure  30.  Sensitivity  of  net  heat  flux  out  and  heat  provided  to  the  ice  from  the  ocean 

C.  VARYING  NET  HEAT  FLUX  WITH  THE  SAME  WIND  SPEED 
1.  Simulation  for  Fifty  Days 

The  effect  of  subtle  changes  in  the  net  heat  flux  was  studied  by  running  the  MICE 
model  for  fifty-days  with  net  heat  fluxes  that  varied  +/-  0.0002  Cal  cm'2  s'1  from  the 
original  forcing  with  the  same  wind  forcing  and  the  same  Greenland  Sea  temperature  and 
salinity  profiles.  The  SHEBA  representative  net  heat  flux  (Figure  14)  was  used  as  a 
baseline  control,  case  1.  For  case  2,  the  net  heat  flux  was  decreased  to  represent  a  smaller 
amount  of  surface  cooling  due  the  presence  of  warmer  air  aloft.  Finally,  in  case  3,  the  net 
heat  flux  was  increased  to  represent  additional  surface  cooling  owed  to  evaporative 
cooling  due  to  the  presence  of  nearby  open  water  or  leads  (Figure  3 1). 
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Net  Heat  Flux 
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Figure  3 1 .  Three  time  series  of  net  heat  flux  forcing  during  simulation. 

As  expected  the  results  show  the  greatest  entrainment  in  case  3  where  the  net 
upward  surface  heat  flux  was  the  largest,  i.e.  most  cooling.  The  increased  cooling  in  case 
3  increased  TKE  and  convection,  which  deepened  the  mixed  layer  23.34  m  over  the  fifty- 
day  simulation.  The  mixed  layer  was  4.09  m  deeper  than  case  1  and  8.69  m  deeper  than 
case  2  (Figure  32).  Case  3  also  fonned  the  most  ice,  3.98  cm,  compared  to  the  3.15  cm 
and  2.21  cm  of  ice  fonned  for  case  1  and  case  2  respectively.  As  a  direct  feedback  to  the 
ice  freezing,  salinity  increased  the  most  in  case  3  due  to  increased  brine  rejection.  Total 
TKE  and  vertical  TKE  for  all  cases  appeared  to  fluctuate  similarly,  due  to  the  identical 
wind  forcing,  but  the  numerical  values  of  TKE  have  differences.  The  maximum  cooling 
case  had  slightly  greater  total  TKE  and  vertical  TKE,  whereas  the  least  cooling  case,  case 
2,  had  the  least  total  and  vertical  TKE.  Figure  33  shows  case  3  had  0.013  cm  s'  greater 
total  TKE  than  case  1  and  0.027  cm2  s'2  greater  total  TKE  than  case  2,  and  0.016  cm2  s'2 
greater  vertical  TKE  than  case  1  and  0.033  cm2  s'2  greater  vertical  TKE  than  case  2, 
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attributable  to  the  increase  in  buoyancy  flux  by  ice  freezing.  The  differences  between 
cases  for  heat  storage  in  the  mixed  layer  was  a  not  result  of  the  temperature  flux  into  the 
bottom  of  the  ice  but  due  to  the  differences  in  temperature  flux  at  the  base  of  the  mxed 
layer.  As  discussed  earlier  the  temperature  flux  into  the  bottom  of  the  ice  varies  with  the 

square  root  of  total  TKE  present  in  the  system,  \[e  .  In  this  case  the  total  TKE  between 
cases  was  similar  and  therefore  the  temperature  flux  into  the  bottom  of  the  ice  was 
similar.  As  a  result,  in  case  3  greater  heat  was  stored  in  the  mixed  layer  given  that  there 
was  greater  entrainment  heat  flux  at  the  base  of  the  mixed  layer. 


Mixed  Layer  Depth  vs.  Time 
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Figure  32.  Varying  net  heat  flux  for  the  same  wind  speed. 
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Vertical  TKE  vs.  Time 


Figure  33.  Total  TKE  and  vertical  TKE  for  vary  net  heat  flux. 

2.  Simulation  for  One-Hundred  and  Fifty  Days 

Running  cases  1,  2  and  3  for  a  longer  period  provides  an  indication  of  the  sea-ice- 
air  interactions  which  may  be  responsible  for  ice  melting  and  the  formation  of  polynyas 
and  leads  despite  a  positive  heat  flux  averaged  throughout  the  period. 
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Figure  31  shows  the  monthly  averaged  net  heat  flux  forcing  during  the  one- 
hundred  and  fifty-day  period.  Throughout  the  model  run,  the  net  heat  flux  was  always 
positive  and  the  system  remained  in  a  deepening  regime.  Ice  begins  to  melt  in  all  cases  at 
Julian  Day  424.  After  Julian  Day  424  the  entrainment  in  all  cases  was  strongly  damped 
due  to  the  negative  buoyancy  flux  in  the  mixed  layer  from  the  stratifying  effect  of  the 
fresh  water  input  due  to  ice  melting  (Figure  34-35).  As  shown  in  the  mixed  layer 
temperature  plot  (Figure  34),  the  mixed  layer  continues  to  provide  heat  to  the  underside 
of  the  ice  during  the  damped  entrainment  periods  by  reducing  the  heat  storage  in  the 
mixed  layer  and  fluxing  it  into  the  ice  bottom.  The  feedback  due  to  ice  melting  is 
illustrated  in  Figure  36  and  Figure  37.  Figure  36  shows  the  ice  transitioning  on  short 
time  scales  between  freezing  and  melting  in  response  to  the  atmosphere  forcing  and 
ocean  feedback.  Figure  37  illustrates  the  natural  damping  effect  on  vertical  TKE,  and 
thereby  entrainment,  due  to  the  ice  melting. 


Mixed  Layer  Depth  vs.  Time  Salinity  vs.  Time 


Temp  vs.  Time  Ice  Thickness  vs.  Time 


Figure  34.  Ice  melting  despite  atmospheric  cooling 
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Total  TKE-Monthly  Forcing  vs.  Time 


Vertical  TKE  vs.  Time 
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Figure  35.  Entrainment  is  damped  in  response  to  the  negative  buoyancy  flux  in  the 

mixed  layer  resulting  from  ice  melting. 
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Figure  36.  Ice  transitioning  on  short  time  scales  between  freezing  and  melting  periods 
in  response  to  the  atmosphere  forcing  and  ocean  feedback. 


Vertical  TKE  vs.  Time 


Figure  37.  Damping  effect  on  vertical  TKE  due  to  the  ice  melting. 
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D.  ATMOSPHERIC  FORCING 


To  address  the  relative  effect  that  temporal  resolution  of  the  forcing  data  has  on 
the  model  skill,  a  moving  average  was  applied  to  the  SHEBA  atmospheric  data  to 
produce  four  time-series,  representative  of  hourly,  daily,  weekly  and  monthly  forcing 
(Figure  34-35).  The  Western  Arctic  and  Greenland  Sea  temperature  and  salinity  profiles 
were  simulated  for  fifty  days  forced  with  each  of  the  four  time-series  of  atmospheric  data. 
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Figure  38.  Moving  average  applied  to  wind  forcing 
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Figure  39.  Moving  average  applied  to  net  heat  flux 

The  plots  from  both  cases  show  the  results  over  the  fifty-day  period  to  be  in 
general  agreement  (Figure  36-37).  However,  there  are  differences  on  a  time  scale  shorter 
than  one  to  two  days.  The  differences  can  be  attributed  to  the  fluctuation  in  TKE 
between  the  different  forcing  time-series  data.  Examining  the  plots  of  total  turbulent 
kinetic  energy  for  the  Western  Arctic  and  Greenland  Sea  model  runs  reveals  that  there  are 
negative  energy  values  predicted  in  both  the  hourly  forced  and  daily  forced  cases  (Figure 
38-39).  Negative  energy  is  unphysical.  Here  TKE  is  being  converted  to  potential  energy 
(PE),  and  the  turbulence  in  the  mixed  layer  cannot  reach  the  base  of  the  earlier 
established  mixed  layer.  This  is  a  period  where  the  system  is  attempting  to  switch  from  a 
deepening  regime  to  a  shoaling  regime.  The  Greenland  Sea,  having  a  deeper  mixed 
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layer,  has  larger  negative  values  of  TKE  because  the  amount  of  shoaling,  or  retreating,  is 
greater. 
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Figure  40.  Western  Arctic  MICE  simulation 
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Figure  4 1 .  Greenland  Sea  MICE  simulation 
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Figure  43.  Greenland  Sea  Total  TKE 

1.  Retreating  Regime 

The  model  can  be  used  as  a  diagnostic  tool  to  predict  the  changing  regimes. 
Understanding  that  water  cannot  be  unmixed  or  “unentrained”,  a  retreating  MICE  model 
was  developed  to  solve  the  system  of  equations  under  a  shoaling  regime.  In  the  retreating 
version  of  the  MICE  model,  entrainment  is  set  to  zero,  and  the  mixed  layer  depth  is 

3 

U 

diagnostically  prescribed  equal  to  the  Obukhov  length  scale,  h  =  L=  *  .  The  effect 

b'w\ 

of  this  is  a  mixed  layer  decoupled  from  the  underlying  ocean.  The  ice  will  still  receive 
heat  from  the  mixed  layer  if  there  was  previous  heat  storage  in  the  mixed  layer,  but  under 
the  shoaling  regime  the  atmospheric  forcing  plays  a  dominant  role. 

The  Greenland  Sea  daily  atmospheric  forcing  case  was  run  for  a  seventeen-day 
period  using  the  retreating  version  of  the  MICE  model  during  the  periods  where  negative 
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TKE  was  predicted.  Figure  40  shows  the  negative  energy  values  occurring  during  three 
periods,  Julian  day:  310,  315,  and  319.  The  models  were  run  to  time  indicated  by  the 
dash  blue  line.  Examining  Figure  41  reveals  similar  system  tendencies  on  the  longer  time 
scale  but  on  the  shorter  period  the  importance  of  changing  regimes  are  clear  and 
significant. 

A  good  example  of  the  decoupling  of  the  heat  source  of  the  deep  ocean  is  at  Julian 
day  315.  Without  the  entrainment  zero,  the  system  tries  to  “unentrain”  temperature  and 
salinity,  which  is  not  physically  realistic.  Furthermore,  the  model  predicts  ice  melting. 
Applying  the  same  forcing  to  a  proper  retreating  model,  no  ice  is  melted  and  the  mixed 
layer  temperature  and  salinity  budgets  are  corrected. 


Greenland  Sea  Daily  Forcing  Total  TKE  vs.  Time 


Figure  44.  Greenland  Sea  Daily  Forcing  Total  TKE 
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Figure  45.  Greenland  Sea  Daily  Forcing  MICE  Model  vs.  Model  verses  Retreating 

Version 
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V.  CONCLUSIONS 


A.  SUMMARY 

A  turbulent  kinetic  energy  (TKE)  budget  model  of  a  coupled  sea-ice-air  mixed 
layer  system  including  realistic  thermodynamics  is  vital  to  the  accurate  prediction  of  ice- 
open  water  boundaries,  polynyas,  and  coastal  and  deep-water  formation  in  the  polar  sea. 
This  work  illustrates  the  complexity  of  sea-ice-air  interactions,  and  the  short  time  scale 
upon  which  the  mixed  layer  responds  to  changes  in  forcing  conditions.  There  is  no  single 
factor  that  is  most  important  in  predicting  high  latitude  sea-ice-air  interactions  because 
the  dynamics  within  the  mixed  layer  are  a  result  of  the  sea-ice-air  system  interaction  as  a 
whole.  The  state  of  the  ice  cover  changes  in  response  to  forcing  from  a  closely  coupled 
atmosphere  and  ocean. 

While  the  atmosphere  generally  governs  the  regime  type,  deepening  or  shoaling, 
the  ocean’s  role  cannot  be  considered  separately.  In  the  deepening  regime  the  ocean 
plays  a  large  role  due  to  entrainment  but  also  plays  a  role  during  shoaling  periods  due  to 
the  heat  storage  capacity  of  the  mixed  layer.  The  MICE  model  allows  a  deeper 
understanding  of  the  limits  and  the  relative  importance  of  the  different  contributors  to  the 
heat,  momentum,  salinity,  and  TKE  budgets.  Prediction  of  buoyancy  flux,  shear 
production,  dissipation,  and  entraimnent  as  functions  of  local  surface  forcing  shed  light 
on  the  intrinsic  feedbacks  within  the  coupled  system. 

It  was  shown  that  changes  in  initial  conditions  such  as  mixed  layer  depth,  and 
temperature  and  salinity  profiles  of  the  polar  oceans  have  significant  effects  on  the 
prediction  of  ice  formation  or  melting.  Increases  or  decreases  in  the  atmospheric  forcing 
may  have  unintuitive  effects  that  only  a  realistic  coupled  system  can  predict.  When 
modeling  ice  thickness  and  the  temperature  and  salinity  profiles  on  short  time  scales,  as 
used  operationally  in  tactically  relevant  problems  such  as  undersea  warfare  or  surface 
warfare,  the  atmospheric  forcing  time  scale  is  crucial  to  accurate  prediction  of  short-term 
regime  changes  that  would  be  missed  by  coarser  resolution  forcing. 

This  work  demonstrates  that  ice  extent  forecast  models  such  as  the  U.S.  Navy 
PIPS  model  potentially  could  be  significantly  improved  by  incorporating  small  time  scale 
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atmospheric  forcing  and  underlying  ocean  structure  in  concert  with  parameterization 
from  the  more  complete  physics  contained  in  the  MICE  model. 


B.  RECOMMENDATIONS  FOR  FURTHER  STUDY 

Currently  the  Matlab®  version  of  the  MICE  model  cannot  automatically 
transition  between  deepening  and  shoaling  regimes  without  user  intervention.  Enhancing 
the  Matlab®  MICE  model  to  diagnose  regime  changes  and  adjust  mixed  layer  depth, 
temperature  profile  and  salinity  profile  would  increase  the  model  capabilities  and  ease  of 
use. 

Further  insight  into  the  nonlinear  physics  and  the  complex  relationships  of  the 
mixed  layer  dynamics  can  be  gained  by  running  the  MICE  model  against  other 
atmospheric  forcing  and  ocean  profile  observations  from  the  Arctic  and  its  marginal  seas. 

The  results  of  this  testing  and  analysis  can  provide  further  guidance  and 
verification  for  the  development  of  the  parameterization  of  heat  flux,  momentum  flux, 
salinity  flux,  and  TKE  budgets  from  a  one  dimensional  model  to  three  dimensional  ocean 
models  to  predict  basin  and  larger-scale  sea-ice-air  interactions. 
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